Molecular & Cellular Proteomics

1Single-cell proteomics (SCP) enables direct measurement of cellular heterogeneity during dynamic biological processes. Here, we applied an SCP workflow to investigate proteome diversity during nerve growth factor (NGF)-induced differentiation of PC12 cells. Differentiated PC12 cells are highly adherent and prone to aggregation, complicating single-cell sample preparation. To address this challenge, sample handling was optimized using gentle dissociation, anti-adhesive conditions, and rapid processing immediately prior to cell isolation. Individual cells were deposited using a refined thermal inkjet (TIJ) dispensing system, enabling accurate single-cell placement with minimal sample loss. Inclusion of the mild nonionic surfactant n-dodecyl-{beta}-D-maltoside (DDM) improved recovery of membrane-associated and other low-solubility proteins. Coupled with high-sensitivity liquid chromatography-ion mobility-mass spectrometry, this workflow consistently quantified approximately 2,000-3,000 proteins per cell across differentiation stages. Single-cell proteomic profiles acquired over the differentiation time course revealed clear separation between undifferentiated and NGF-treated cells by Day 6. At later stages (Days 4-6), cells further partitioned into two distinct subpopulations with protein expression patterns not evident in bulk measurements. Dimensionality reduction and non-negative matrix factorization identified multiple proteomic states coexisting within the same differentiation stages, characterized by coordinated differences in pathways related to intracellular trafficking, protein translation, and neuronal structural organization. Together, these results show that while global proteome remodeling during PC12 differentiation is captured in both bulk and single-cell data, single-cell proteomics uniquely resolves functionally distinct cellular subpopulations that are masked in population-averaged analyses.

The unfolded protein response (UPR) helps reinstate cellular proteostasis upon an accumulation of misfolded proteins in the endoplasmic reticulum (ER), in part through ER-associated degradation (ERAD). Ube2j2 is an ER-localized E2 ubiquitin-conjugating enzyme that participates in ERAD. We used mass spectrometry analysis of cultured U2OS cells to investigate how the loss of Ube2j2 affects the cellular proteome in response to tunicamycin-induced ER stress. We constructed a network of twelve statistically distinct modules of protein abundance profiles across conditions. We describe the Gene Ontology annotations for each module along with the "hub gene" proteins whose abundance levels most closely adhere to each modules protein abundance profile. Our analysis identifies known Ube2j2-associated pathways (e.g., the UPR and ERAD) and cellular functions that were previously unassociated with Ube2j2 (e.g., RNA metabolism, ER-Golgi transport, and cell-cycle progression). These data are available via ProteomeXchange with identifier PXD076153 and provide avenues for further investigation into the cellular functions of Ube2j2 under basal and ER-stressed conditions.

Understanding allelic variability is crucial for elucidating intrinsic bacterial mechanisms and distinguishing phenotypic profiles. However, such variability poses a major challenge for the reliable identification of proteins in data-independent acquisition (DIA) proteomics. To address this, we developed an analytical workflow that integrates protein sequence variability to enhance proteome coverage. Fifteen Legionella pneumophila isolates were analyzed using DIA-NN, with spectral libraries generated either from a reference proteome or incorporating allelic variability. Our workflow includes protein clustering and subsequent protein inference from these clusters, allowing the accurate assignment of shared and variant-specific peptides. Integration of variability enabled the identification of a comparable number of proteins as the reference proteome while capturing between 28 and 77 % of variant-specific sequences in each isolate, all while maintaining a low false positive rate. These findings demonstrate that accounting for allelic variability substantially improves proteomic coverage and identification confidence, providing a more comprehensive view of the proteome. This approach facilitates a deeper understanding of biological mechanisms and enables precise bacterial proteotyping of Legionella pneumophila isolates.

A recent clinical study tested the effects of two different monoclonal antibodies (mAbs) (siltuximab, anti-IL6; tocilizumab, anti-IL6R) on the fate and function of T-cells in people with type 1 diabetes. While both mAbs affect the response of T-cells to stimulation, they have very different, sometimes opposing mechanisms. Here, we use mass-spectrometry based proteomics to analyze longitudinal serum samples (baseline and two weeks post-treatment) from 20 clinical trial participants to examine the effects of siltuximab and tocilizumab on extracellular vesicles. To accomplish this, serum samples were enriched for extracellular vesicles with Mag-Net and analyzed by LC-MS/MS to identify significantly differentially abundant protein groups and pathways. Proteome analysis confirmed highly reproducible measurements across multiple draw dates. In total, we quantified >3300 protein groups of which 46 protein groups had significantly altered abundance after mAb treatment. Tocilizumab altered pathways associated with proteostasis (neddylation) and pre-notch transcription and translation. Siltuximab altered FCGR activation pathway members. In addition, quantitation of the monoclonal antibody therapies themselves enabled the measurement of the correlation between drug amounts and impacted proteins. Taken together, this work demonstrates the utility of the Mag-Net method to evaluate the impacts of therapeutic interventions on serum extracellular vesicles.

Combination therapies are widely used in acute myeloid leukemia (AML), but systematic datasets capturing proteome-wide responses to multi-drug perturbations remain limited. Here we present CoPISA (Combinatorial Proteome Integral Solubility/Stability Alteration), a quantitative proteomics assay designed to profile protein solubility and stability responses to single and combined drug treatments. The dataset includes two AML drug pairs (LY3009120-sapanisertib and ruxolitinib-ulixertinib) applied to four AML cell lines (MOLM-13, MOLM-16, SKM-1, and NOMO-1) under control, single-agent, and combination conditions in both lysate and intact-cell formats. Thermal solubility profiling coupled with TMT-based multiplexed LC-MS/MS generated 16 TMT16-plex experiments comprising 192 LC-MS/MS raw files, providing deep proteome coverage across treatments and biological contexts. The resource includes raw and processed proteomics data, detailed experimental metadata in Sample and Data Relationship Format (SDRF), and reproducible analysis scripts for reporter normalization, protein-level aggregation, statistical modeling, and classification of combinatorial response patterns. The experimental design enables identification of proteins responding uniquely to combination treatments as well as overlapping single-agent effects. Technical validation demonstrates reproducible quantification across multiplex experiments and assay formats. All data are publicly available through the PRIDE repository (PXD066812) together with analysis code, enabling independent reanalysis and method development. This dataset provides a benchmark resource for studying proteome responses to drug combinations, comparing lysate and intact-cell perturbation profiles, developing computational approaches for combinatorial target inference, and supporting training in computational proteomics.

Large-scale proteomics relies heavily on target-decoy competition for false discovery rate estimation in peptide identification, and the performance of this strategy depends strongly on the design of the decoy database. Classical generators such as reversal and shuffling remain widely used. Here, we introduce protein language model-based (PLM) decoy generation for peptide identification and benchmark it against classical strategies. We evaluate these approaches using three complementary quality-control layers: sequence-based separability, search-engine-agnostic spectral-space diagnostics, and end-to-end mass spectrometry benchmarks, including pipelines with rescoring. Across these analyses, PLM-based decoys are harder for sequence-only neural networks to distinguish than most classical generators, suggesting fewer obvious sequence-level artifacts. However, this signal is only weakly informative for search performance. Spectral diagnostics further show that short peptides occupy a particularly crowded target-decoy space and are therefore especially prone to local collisions across all generators. In full search pipelines, reverse decoys remain a strong baseline, and current PLM-based generators do not yet provide a clear overall advantage. We therefore view PLM-based decoys not as universal replacements for reverse decoys, but as tunable tools for benchmarking, diagnostics, stress testing, and future adaptive decoy optimization, with increasing value as search models become more expressive.

Single-cell proteomics (SCP) by mass spectrometry can now quantify hundreds to thousands of proteins per cell, but the field still lacks standardised analytical pipelines that accommodate the diversity of instruments, sample preparation workflows and biological contexts encountered in practice. Existing workflows, largely adapted from single-cell transcriptomics, do not account for the informative missingness, pervasive ambient protein contamination and limited feature space that distinguish proteomic from transcriptomic data. In addition, cell type annotation remains a manual bottleneck that is subjective, difficult to reproduce and hard to scale. Here we present an end-to-end pipeline that integrates adaptive quality control, entropy-guided iterative batch correction, multi-modal marker discovery that exploits detection patterns unique to proteomics, and context-aware annotation by large language models (LLMs) coupled to structured contradiction reasoning and orthogonal data-driven validation. Benchmarking on published single-cell proteomic datasets from developing human brain and glioblastoma-associated neutrophils revealed systematic LLM failure modes, including context-insensitive marker vocabulary and misinterpretation of phagocytic or lytic cell states. We addressed these errors using a three-round prompt architecture that combines general biological principles with auto-generated dataset-specific constraints. In held-out validation on a skin tumour dataset acquired, the pipeline showed high concordance with FACS-sorted ground truth. In the caerulein-injured pancreas, orthogonal immunohistochemistry further supported annotations of macrophage, stellate and immune populations. The pipeline is fully automated under fixed settings, and available as Context-Aware Single-Cell Proteomics Analysis (CASPA), providing SCP laboratories and facilities with a reproducible workflow that delivers interpretable, confidence-quantified annotations suitable for downstream expert review.

Peptide-level analyses are becoming increasingly popular in mass spectrometry-based proteomics and are being applied, for example, in immunopeptidomics, structural proteomics, and analyses of post-translational modifications. In such analyses, peptides that are not biologically meaningful but instead arise as artifacts prior to mass spectrometry analysis pose the risk of data misinterpretation. Here, we describe an approach based on retention time analysis and precise chromatographic peak matching to identify peptides generated by in-source fragmentation (ISF), which occurs between chromatographic separation of peptide mixtures and the first mass filter of a tandem mass spectrometer (MS). To understand the prevalence and properties of ISF, we generated 13 proteomics datasets and analyzed them along with additional 25 previously published datasets spanning a broad range of sample types, MS, and proteomics approaches including classical bottom-up proteomics, immunopeptidomics, structural proteomics, and phosphoproteomics. We found that, in typical trypsin-digested samples on average 1 % of fully-tryptic peptides and 22 % of semi-tryptic peptides originated from ISF. However, we observed large variations between datasets, and in-source fragments exceeded, in some cases, a third of the total peptide identifications. The extent of ISF was dependent on the peptide sequence, the instrument, method parameters, and sample complexity. Although ISF did not impair relative quantification across samples, it generated peptides that could be misinterpreted qualitatively, inflated peptide identifications, and comprised up to 37 percent of peptides shorter than 9 amino acids in immunopeptidomics datasets. We propose that, for peptide-centric applications, our open-source ISF detection approach be used to re-annotate peptides generated by ISF and remove them to avoid misinterpretation of data. ISF is an increasing concern with improving mass spectrometers, as they enable detection of an ever-increasing number of m/z features, including low abundance features like ISF products. Our work thus addresses a growing issue in proteomics and presents solutions to mitigate the impact of in-source fragment peptides. In the future, improved feature detection algorithms may enable elucidation of new ISF patterns affecting side chains that have been missed so far, which could contribute to explaining the vast space of as-yet unannotated proteomics data.

Open modification search (OMS) strategies have gained popularity in mass spectrometry-based proteomics for identification of peptides carrying unknown or unexpected post-translational modifications. However, most OMS search engines report only the overall mass difference between the precursor and the matched peptide and do not explicitly identify or score combinations of multiple modifications at the peptide-spectrum match (PSM) level, leaving the interpretation of mass shifts up to the end user and to using downstream analysis tools. Here, we introduce MS Andrea, a novel OMS search engine developed to directly identify and score combinations of up to four variable modifications per peptide without having to predefine them. MS Andrea uses a sequence tag-based strategy to efficiently filter candidate peptides prior to scoring. Remaining candidates are evaluated using the MS Amanda scoring function, first considering fixed modifications only, followed by a second scoring stage in which combinations of modifications from the Unimod database are considered based on the observed mass difference and matched to the spectrum. We evaluated MS Andrea using phosphopeptide datasets from HeLa cells and Arabidopsis thaliana and compared its performance with the widely used OMS engines MSFragger and Sage. Across datasets, MS Andrea identified the highest number of PSMs at 1% false discovery rate while achieving comparable peptide-level identifications. Importantly, MS Andrea directly reports modification identities and sites at the PSM level and enables the identification of peptides having up to four variable modifications. Together, these results demonstrate that MS Andrea facilitates more detailed and interpretable characterization of peptide modifications while maintaining competitive identification performance in OMS-based proteomic analyses. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=132 SRC="FIGDIR/small/714851v1_ufig1.gif" ALT="Figure 1"> View larger version (19K): org.highwire.dtl.DTLVardef@52f65forg.highwire.dtl.DTLVardef@acf4e3org.highwire.dtl.DTLVardef@10171caorg.highwire.dtl.DTLVardef@1d594ad_HPS_FORMAT_FIGEXP M_FIG C_FIG

SummaryExtracellular vesicles (EVs) are bilayer vesicles that carry a diverse cargo of molecules, such as nucleic acids, proteins and metabolites. These EVs can be transported throughout the organism to specific recipient tissues. For this reason, EVs have been recognized as pivotal mediators of cell-to-cell communication (CCC). Importantly, alterations in EV-mediated communication have been linked to pathological processes, further highlighting their biological relevance. However, the in silico exploration of the functional effects of EV cargo in recipient tissues remains limited due to the lack of dedicated tools that can be applied to EV omics datasets. Most current bioinformatics tools for assessing CCC rely on ligand-mediated communication and therefore cannot be used to explore EV-mediated communication. To address this gap, we developed EV-Net, a bioinformatics tool designed to explore the effects of EV cargo on recipient tissues. EV-Net was built by adapting NicheNet, a CCC bioinformatics tool that relies on ligand-receptor mediated communication, for the analysis of EVs proteomics and RNA-seq data. The EV-Net framework enables the identification and prioritization of EV cargo molecules with high regulatory potential in a recipient tissue of interest. This prioritization facilitates the systematic translation of EV cargo profiles into testable biological hypotheses. Availability and documentationThe source code of EV-Net is stored in GitHub https://github.com/torrejoNia/EV-Net alongside instructions on how to install it. Comprehensive tutorials and additional documentation are available at https://torrejonia.github.io/EV-Net/. The datasets used in the use cases were deposited in Zenodo. The corresponding Zenodo links are provided in the tutorials for each use case. This software is distributed under a GLP3 licence.

Astrocytes are a subset of glial cells in the central nervous system (CNS) that support numerous processes essential for brain function. Their functional diversity is thought to arise from specialized subpopulations with distinct molecular profiles. Although single-cell and single-nucleus RNA sequencing (scRNA-seq and snRNA-seq) have greatly advanced our understanding of astrocyte transcriptomic heterogeneity, mRNA abundance does not always correlate with protein levels because of post-transcriptional and translational regulation. Therefore, studying protein profiles remains essential to accurately capture astrocyte functional states and heterogeneity. Here, we used Microscoop Mint, a microscopy-guided spatial proteomics platform that integrates subcellular, region-specific sample preparation with LC-MS/MS-based mass spectrometry, enabling direct protein profiling of astrocytes in paraformaldehyde-fixed, optimal cutting temperature (OCT)-embedded mouse brain tissue. By applying this approach, we uncovered distinct regional-associated astrocyte proteomic signatures in the cerebral cortex and hippocampus and selected novel candidate protein markers for subsequent validation by immunofluorescence. Notably, MINK1 and PLEKHB1 showed preferential expression in hippocampal and cortical astrocytes, respectively, highlighting their potential as region-specific astrocyte markers. Overall, this strategy enables high-precision, unbiased spatial proteomic discovery at subcellular resolution, providing a powerful framework for linking molecular diversity to functional specialization in astrocyte biology.

Cells contain different organelles and compartments that are essential for cellular function and life. These organelles and compartments need to communicate to assess cellular state in a changing environment, adapt to the new situation, and also to ensure functionality and homeostasis. Moreover, organization and communication differ between cell types. However, our knowledge about these changes is still rather scarce. Subcellular spatial proteomics aims to fill this knowledge gap. While proximity labeling techniques represent a great advance, they do not provide precise spatial resolution. To overcome this limitation, we developed SPEx (Subcellular spatial Proteomics coupled to Expansion), in which we first expand cells about 10- fold, laser micro-dissect regions of interests and then perform mass spectrometry-based proteomics on these samples. We demonstrate the effectiveness of SPEx by determining the proteome of the Golgi, the nucleus and nucleoli. Satisfyingly, we also identify novel components of these organelles. Combining inexpensive already existing technologies makes SPEx readily usable by the wider scientific community.

Mass spectrometry imaging (MSI) characterizes the spatial heterogeneity of molecular abundances in biological samples. Experiments with complex designs, involving multiple conditions and multiple samples, provide particularly useful insight into differential abundance of analytes. However, analyses of these experiments require attention to details such as signal processing, selection of regions of interest, and statistical methodology. This manuscript contributes a statistical analysis workflow for detecting differentially abundant analytes in MSI experiments with complex designs. Using a case study of histologic samples of human tibial plateaus from knees of osteoarthritis patients and cadaveric controls, as well as simulated datasets, we illustrate the impact of the analysis decisions. We illustrate the importance of signal processing and feature aggregation for preserving biological relevance and alleviating the stringency of multiple testing. We further demonstrate the importance of selecting regions of interest in ways that are compatible with differential analysis. Finally, we contrast several common statistical models for differential analysis, showcase the appropriate use of replication, and demonstrate model-based calculation of sample size for followup investigations. The discussion is accompanied by detailed recommendations and an open-source R-based implementation that can be followed by other investigations.

Data-independent acquisition (DIA) proteomics enables reproducible and systematic peptide detection and quantification, and trapped ion mobility spectrometry (TIMS) on the timsTOF platform further improves DIA by synchronizing ion mobility separation with quadrupole precursor sampling. Analyzing the highly multiplexed spectra generated by DIA typically relies on spectral libraries, and fully leveraging the additional ion mobility dimension requires these libraries to include accurate retention time, fragment ion intensity, and ion mobility annotations. Existing in silico spectral library generation tools either lack ion mobility support entirely or rely on models trained on data-dependent acquisition (DDA) data, that can introduce a mismatch that may not capture unique experiment-specific biases when applied to each respective timsTOF dataset. Carafe is a software tool that uses deep learning models to generate high-quality, experiment-specific in silico libraries by training directly on DIA data. In this study, we extend Carafe to generate libraries for timsTOF DIA data, which involves fine-tuning retention time (RT), fragment ion intensity, and ion mobility prediction models using timsTOF DIA data. Carafe2 operates directly on native timsTOF raw data (Bruker .d directories) without the need for data conversion. We demonstrate the performance of Carafe2 across a wide range of DIA applications, including global proteome, phosphoproteome, and plasma proteome datasets. Comparing Carafe2 fine-tuned RT, fragment ion intensity, and ion mobility prediction models with pretrained DDA models, we find that Carafe2 models outperform pretrained models on a variety of DIA datasets. We then demonstrate the utility of in silico libraries generated by Carafe2 for peptide detection on several different types of timsTOF DIA datasets by comparing with the libraries generated with DDA-trained AlphaPeptDeep models, DIA-NN built-in models, and empirical spectral libraries generated from DDA experiments.

Limited proteolysis coupled to mass spectrometry (LiP-MS) is a peptide-centric conformational proteomics approach during which a brief incubation with a non-specific protease (e.g., proteinase K) under native conditions generates structural fingerprints that report on treatment-induced conformational changes, which is followed by a tryptic digest under denaturing conditions allowing to read out these fingerprints 1. In contrast, the recently introduced peptide-centric local stability assay (PELSA) uses a high trypsin-to-substrate ratio under native conditions to release fully tryptic peptides that reflect structural stability upon ligand binding 2. In their paper, Li et al. compared PELSA and LiP-MS across several benchmarks and reported that PELSA exhibited quantitative sensitivity comparable to or exceeding LiP-MS. Notably, PELSA quantified a 21-fold greater rapamycin-induced change for FKBP1A compared to LiP-MS. Because such claims influence method selection for conformational proteomics, we reanalyzed the publicly deposited datasets underlying these comparisons and assessed the experimental and analytical choices that contributed to the reported effect sizes. Our evaluation indicates that the reported 21-fold difference arises from non-matched experimental conditions and undisclosed data imputation, and that conclusions regarding quantitative superiority or biological interpretability should therefore be treated with caution.

SummaryProteoPy is a lightweight Python library for protein- and peptide-level quantitative proteomics analysis, built around the AnnData class as its core data structure. It streamlines data import, preprocessing, and differential analysis while preserving all metadata within a single object. A reimplementation of our previously published COPF algorithm enables proteoform group inference directly from peptide-level data, facilitating the identification of proteoform-specific regulation and isoform usage. Designed for accessibility and flexibility, ProteoPy simplifies analysis for non-specialists and provides an extensible foundation for advanced proteomics workflows, seamlessly integrating with the scanpy and muon ecosystems for reproducible and scalable multi-omics analysis. Availability and implementationProteoPy is implemented in Python 3 and publicly available on GitHub: https://github.com/UKHD-NP/proteopy under the Apache 2.0 license. Contactisabell.bludau@med.uni-heidelberg.de Supplementary informationTutorial notebooks for ProteoPy are included as supplementary data and are also available on GitHub: https://github.com/UKHD-NP/proteopy/tree/main/docs/tutorials.

Defining how proteins change over developmental time is amenable to studies deciphering regulatory genetic networks in vertebrate development, biology, and pharmacology. In an approach toward such quantitative studies of dynamic network behavior, we produced an atlas using the mass spectrometry-based method to investigate protein expression changes across 16 time points from the zygote to the early pharyngula stage zebrafish embryos. We systematically summarize 8 clusters for interrogating changes in protein expression associated with the development of zebrafish embryos. Specifically, we identified a class of zinc finger-related transcription factors primarily located on the long arm of chromosome 4, which are highly expressed during zygotic genome activation. Furthermore, we highlight the power of this analysis to assign developmental stage-specific expression information to chromosomes and tissues. Time-resolved analyses reveal significant discordance between differential transcript and protein expression, whereas no time lag is observed for proteins involved in stable and fundamental biological processes, such as metabolism (e.g., Ppt2a and Gatm), cytoskeletal organization (e.g., Col18a1), and the translation machinery (e.g., Eif4enif1). This atlas offers high-resolution and in-depth molecular insights into zebrafish development, providing a resource for developmental biologists to generate hypotheses for functional analysis of proteins during early vertebrate embryogenesis. HighlightsO_LIA global protein expression database with high time resolution is created for zebrafish embryos. C_LIO_LIDistinct patterns of protein expression correlate with biological processes. C_LIO_LITranscription factors have a burst of expression from the gastrulation stage. C_LIO_LIDevelopmental stage-specific protein expressions were assigned to chromosomes and tissues. C_LIO_LIHigh-resolution embryonic transcriptome and proteome datasets were compared and connected. C_LI

Effective predictive biomarkers for immune checkpoint inhibitor (ICI) therapy remain an unmet need across solid tumors. Here, we present an integrated spatial proteomics workflow that combines in situ proximity ligation assay with multiplexed immunofluorescence to directly resolve PD1/PDL1 signaling events at the level of defined cellular phenotypes and their spatial organization within intact tumor tissue. Applied as a proof of concept to tumor samples from patients with metastatic urothelial carcinoma treated with pembrolizumab, this approach reveals that PD1/PDL1 interactions specifically involving cytotoxic CD8CD3 T cells are significantly enriched in complete responders, while such interactions are rare in patients with progressive disease. This interaction defined T cell subset achieves superior discrimination of clinical response compared to single marker PDL1 expression or immune cell abundance alone. By integrating direct detection of protein protein interactions with high dimensional single cell phenotyping, our workflow provides a mechanistically informed, spatially resolved biomarker of functional immune engagement. Beyond urothelial carcinoma, this platform establishes a generalizable framework for translating spatial signaling biology into predictive tools for immunotherapy response across tumor types.

Accurate genome annotation is fundamental to modern biology, yet distinguishing authentic protein-coding sequences from prediction artifacts remains challenging, particularly in complex plant genomes where automated methods are error-prone and manual curation is rarely feasible due to prohibitive time and costs. Here, we present GAP-MS (Gene model Assessment using Peptides from Mass Spectrometry), an automated proteogenomic pipeline that leverages mass spectrometry evidence to systematically validate the protein-level accuracy of predicted gene models. Applied across 9 major crop species, GAP-MS consistently improved prediction precision for four widely used gene prediction tools. In addition to filtering erroneous models, the pipeline identified hundreds of previously missing gene models from current standard reference annotations. These peptide-supported loci were further verified by transcriptional evidence, well-supported functional annotations, and high coding-potential scores. Together, these results demonstrate that direct proteomic evidence provides a robust framework for resolving annotation ambiguities, defining high-confidence reference proteomes, and uncovering overlooked protein-coding genes, while facilitating the identification of sequences that may require further investigation. GAP-MS is freely available as a web interface at https://webclu.bio.wzw.tum.de/gapms/.

Lymph nodes (LN) are key secondary lymphoid organs (SLO) for a coordinated immune response. They have been extensively characterized by numerous investigative techniques chiefly as single cell suspensions because they are composed of vagile yet crowded hematolymphoid elements, unfriendly to spatial tissue organization-saving techniques. We comprehensively classify in situ all cells of 19 human LN free of pathology with a 78-marker antibody panel, an hyperplexed cyclic staining method, MILAN, and an analytical bioinformatic pipeline, BRAQUE. A total of 77 cell types were classified, encompassing T, B, innate immune and stromal cells. CD4 and CD8 T-cells were classified into 27 unique subsets by leveraging the expression profiles of TCF7, the presence of co-inhibitory receptors and the spatial distribution. CD5 and TCF7 expression defined novel B-cell types. CD27+ mature B-cells occupied previously unrecognized nodal spaces non-overlapping with the cortex and the plasma-cell rich medullary cords. Type 2 conventional dendritic cells were located in nodular paracortical aggregates. Statistically controlled pairwise neighborhood analysis showed sparse cell-cell interactions, known and new neighbors, established and novel LN landscape niches. A high-dimensional proteomic interrogation of the normal human LN provides spatial allocation of known cell types, novel interactions and the landscape organization.